Screening and Dissipation at the Superconductor-Insulator Transition 
Induced by a Metallic Ground Plane 
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We study localization phenomena in two dimensional systems of charged particles in the presence 
of a metallic ground plane with a particular focus on the superconductor-insulator transition. The 
ground plane introduces a screening of the long-range Coulomb interaction, and provides a source 
of dissipation due to the gap less diffusive electrons. The interplay of these two effects leads to 
interesting physical phenomena which are analysed in detail in this paper. We argue that the generic 
superconductor-insulator transition of charged particles in the presence of the ground plane may be 
controlled by a. fixed line with variable critical exponents. This is illustrated by an explicit calculation 
in an appropriate e expansion. In contrast, the universal properties of the superconductor-Mott 
insulator transition in the clean limit at commensurate densities are shown to be unmodified by 
either the metal or the long-range Coulomb interaction. A similar fixed line can arise in the presence 
of a metallic ground plane for quantum Hall plateau transitions. Implications for experiments on 
Josephson-junction arrays and quantum Hall systems are described. 



I. INTRODUCTION 

Despite many decades of effort, the properties of mat- 
ter in the vicinity of various transitions from conducting 
to insulating phases remains poorly undcrstod. The best 
known example is the transition from a metal to an in- 
sulator. During the last few years, a number of such 
phenomena have been subjected to serious experimental 
study, particularly in low dimensional systems. These in- 
clude transitions from a superconductor to an insulator 
[1] , various transitions in two dimensional systems show- 
ing the quantum Hall effect [2], and an apparent tran- 
sition from an as yet unidentified conducting phase to 
an insulator in two dimensional Silicon MOSFET's [3]. 
This increase in experimental activity has raised ques- 
tions that require an improved theory of such localization 
phenomena. 

In this paper, we study localization phenomena in 
two dimensional systems in the presence of a proximate 
ground plane formed by a diffusive two dimensional elec- 
tron gas. Our focus will be the superconductor-insulator 
transition, though many of our results generalize to other 
transitions as well. Recent experiments have studied the 
effect of such a metallic plane on the quantum transition 
to superconductivity in Josephson junction arrays [4] and 
in homogcnously disordered thin films in a magnetic field 
[5] . In these experiments, the actual tunneling of charged 
excitations between the superconductor and the metallic 
plane is very weak, and may be ignored. However, the 
charge carriers in the two systems are coupled together 
by the Coulomb interaction. Such a metallic plane has 
several effects. First, one might expect some screening of 
the long-ranged Coulomb interaction between the charges 
in the system undergoing the localization transition. Sec- 
ond, the coupling to the gapless excitations in the metal 



provides a source of "dissipation" . The interplay of these 
two effects leads to interesting differences in the localiza- 
tion phenomena from that in the absence of the metallic 
plane. In particular, we argue that the zero temperature 
localization transition may be controlled by a fixed line 
with variable critical exponents. 

The interest in studying the role of a metallic plane 
stems from several different directions. A particularly 
difficult issue for localization theory is the role of the 
long-range Coulomb interaction and its possible screen- 
ing. For the integer quantum Hall transitions, for 
instance, while short-range interactions have been ar- 
gued to be irrelevant at the non-interacting fixed point 
[6], the nature of the criticality with Coulomb interac- 
tion is very poorly understood. For metal-insulator or 
superconductor-insulator transitions, the Coulomb inter- 
action is screened in the conducting phase but retains 
its long-range character in the insulating phase. This 
evolution in the screening properties of the system as 
it moves through the transition contributes to the diffi- 
culty in developing a theoretical understanding. Intro- 
ducing a metallic plane enables exerting some control on 
the screening and may enable separating out phenomena 
that are specific to the long-range Coulomb interaction. 

A different motivation comes from theoretical attempts 
to understand experiments on the destruction of super- 
conductivity at zero temperature due to quantum ef- 
fects. In situations where the transition is to an insu- 
lator, and at asymptotically low temperatures, a "boson- 
only" model that describes the localization of Cooper 
pairs while ignoring the fermionic degrees of freedom 
is expected to describe the universal critical properties 
[7]. However, in a number of experiments, the situa- 
tion is different [1,5]. The non-superconducting phase 
may be a weakly localized metal, or in some cases ap- 
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pear to be truly metallic (with no trace of localization 
effects) in the temperature range currently probed in ex- 
periments. In these cases, a "boson-only" model is pre- 
sumably inadequate, and it is ncccsary to include the 
underlying electrons. However a full treatment of the 
problem of Cooper pairs coupled to gapless electrons is 
prohibitively difficult it therefore helps to study simpler 
situations where some if not all the effects are treated re- 
liably. Such a simplification is provided by considering 
"boson-only" models coupled to gapless metallic elec- 
trons through density-density Coulomb interactions, in 
which the decay of the boson into pairs of electrons is ig- 
nored. Precisely this situation is realized in experiments 
on Josephson junction arrays where a proximate metallic 
plane is introduced [4]. Remarkably, as we show in this 
paper, such a coupling between the Cooper pairs and the 
metallic electrons has a profound effect on the properties 
of the transition. 

The nature of the superconducting phase and its tran- 
sition to the insulator in the presence of a metallic plane 
has been addressed before in the literature. We believe 
that these previous treatments miss several crucial as- 
pects of the physics. Ref. [8] considered a model of 
a Josephson junction array with a short-ranged capaci- 
tance matrix that essentially amounts to assuming a loga- 
rithmic Coulomb interaction energy between the charges 
rather than the 1 /r interaction that obtains for long dis- 
tances and controls the universal properties of the tran- 
sition. This difference profoundly modifies the physics — 
hence the results of Ref. [8] are not capable of describ- 
ing the experiments of interest [4,5]. There arc also a 
number of papers on "local dissipative" models, where 
the dissipation is claimed to arise from spatially localized 
excitations. In a short section below, we discuss the rela- 
tionship between these models and our work. First in our 
work unlike in Ref. [9] , the physical origin of the dissipa- 
tion is clear - it is due to the gapless diffusive electrons in 
the metal. Furthermore, as we disciiss in detail, the fixed 
line obtained here has an entirely different origin from the 
fixed line claimed to exist in the model studied in Ref. [9] 
(see also [10]). Deep inside the superconducting phase, 
the dynamics of the phase of the superconducting order 
parameter in the presence of the coupling to the metal 
was considered by Gaitonde [11], who argued that the 
plasmon was overdamped at long wavelengths. We de- 
rive the correct phase dynamics and show that in fact the 
plasmon mode survives at long wavelengths. 

For most of this paper, we focus on the 
superconductor-insulator transition of repulsively inter- 
acting bosons in the presence of density-density coupling 
to a diffusive two dimensional electron gas. We begin 
by first considering the effects of the metallic plane on 
the superconducting phase. We show that the plasmon 
mode of the superconductor survives essentially unmod- 
ified. This is contrary to naive expectations based on 
the assumption that a short-ranged screened interaction 



between the bosons correctly describes all aspects of the 
physics with the metal present. We find that, because 
of their slow diffusive motion, the electrons in the metal 
are unable to screen out the faster charge fluctuations 
associated with the plasmon. This analysis sets the stage 
for a discussion of the general properties of the transi- 
tion to the insulator. With the metallic plane present, 
we give a general argument that points to the existence 
of a fixed line controlling the transition. The fixed line 
is parametrized by the value of the conductivity of the 
metal. Universal critical properties will then be deter- 
mined by the metal conductivity. 

To substantiate our general arguments we present a 
detailed analysis of a specific model which is appropri- 
ate to describe the experiments in Ref. [4]. We consider 
the superconductor-insulator transition at commensurate 
density for bosons on a square lattice in the presence 
of the metal. We first provide crude estimates of the 
phase boundary and show that the presence of the metal 
stabilizes the superconducting phase relative to the in- 
sulator. We derive an effective action for the Cooper 
pair degrees of freedom by integrating out the diffusive 
metallic electrons which is then used to discuss the crit- 
ical phenomena at the transition. In the clean system, 
we argue that neither the metal nor the long-ranged part 
of the Coulomb interaction affects the critical properties 
so that the transition is in the universality class of the 
D = 2 -\- 1 XY model. This result is strongly suggested 
by two separate calculations — a large- iV generalization of 
our effective model, and an e expansion of an appropri- 
ate generalization to dimensions other than 2. However 
moving away from this idealized limit is expected to al- 
ter the critical properties and produce the fixed line ex- 
pected from the general arguments advanced earlier. We 
illustrate this by including disorder which nevertheless 
preserves the particle/hole symmetry of the commensu- 
rate problem. A double e expansion [12] along the lines of 
Ref. [13] is then performed, including both the Coulomb 
interaction and the effect of the metal, which indeed finds 
the expected fixed line. 

The generality of our arguments for the existence of 
a fixed line in the presence of the metal is then illus- 
trated in the context of other concrete yet tractable mod- 
els of localization transitions. Specifically, we show that 
a model for the quantum Hall transition considered in 
Ref. [14] should also possesses this feature. Indeed we 
expect that the generic quantum Hall transition with 
long-range Coulomb interactions could, in the presence 
of a proximate metal, also be controlled by a fixed line 
with variable critical exponents. 

The implications of our results for experiments on 
Josephson-j unction arrays, superconducting films, and 
quantum Hall transitions are considered in Section VII. 
An important qualitative lesson is that the presence of a 
proximate metallic plane does not by itself guarantee that 
the universality class is that of the short-ranged prob- 
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1cm. Under what circumstances then is the latter real- 
ized? For the generic superconductor-insulator or quan- 
tum Hall transition, our analysis indicates that when the 
fixed line scenario is realized, and the conductivity of the 
metal is very large (compared to e^//i) there will be a 
region near the transition where the universal behaviour 
is controlled by the short-range fixed point. Very close 
to the critical point however there will be a crossover 
to behaviour controlled by the appropriate point on the 
fixed line. The presence of a fixed line directly impacts 
experiments seeking to measure universal transport at 
these two dimensional localization transitions. Indeed, a 
temperature-independent conductivity will obtain right 
at the transition point, with a value that will vary along 
the fixed line and hence may appear to be non-universal. 
However, this conductivity will be universally related to 
the critical exponents (which too will vary along the fixed 
line). Experimental realization of the fixed line may be 
easier to observe for the quantum hall transition than for 
the superconductor-insulator transition as the universal 
scaling regime appears to be more easily accesible in the 
former. 

It is interesting to speculate on the implications of 
our results for the general superconductor-metal tran- 
sition. Our analysis shows that a simple Coulomb cou- 
pling between the order parameter (boson) degrees of 
freedom and gapless diffusive electrons already leads to 
profound modifications of the universal properties from 
that of a boson-only model. To correctly describe the 
superconductor-metal transition , it is necessary to in- 
clude processes where the Cooper pairs can decay into a 
pair of electrons and vice versa into the model. If any 
portion of the fixed line survives this inclusion, then the 
conductivity at the transition will again appear to be 
non-universal and will be determined by the value of the 
marginal coupling at the point on the fixed line that con- 
trols the transition. Wo note that experiments on two- 
dimensional superconductors undergoing a transition to a 
weakly localized metal do see a temperature- independent 
but apparently non-universal conductivity at the transi- 
tion. 

II. GENERAL ARGUMENTS 

We begin with some simple observations on the effects 
of a density- density coupling between the system of in- 
terest and a diffusive 2DEG. For concreteness, we con- 
sider the physical situation shown in Figure 1. A two- 
dimensional Josephson junction array ( JJA) is separated 
by an insulating slab from a two dimensional electron 
gas (2DEG), with conductivity Ue- The Cooper pairs in 
the JJA couple to the electrons in the 2DEG solely via 
the Coulomb interaction. We assume that the electron 
motion in the 2DEG is diffusive, and ignore weak local- 
ization effects. For the experimental situation of interest 



this is well justified since the typical 2DEG conductiv- 
ities are large enough that weak localization effects are 
insignificant at currently accessible temperatures. 




Insulator 



FIG. 1. The setup considered in this paper, based on the 

experiments of Rimberg et al. [4] A Josephson junction array 
(JJA) is separated from a two dimensional electron gas (with 
conductivity Ue that can be tuned) by an insulating layer that 
blocks tunneling of electrons. The Cooper pairs in the JJA 
therefore interact with the electrons in the 2DEG solely via 
the Coulomb interaction. 

A. Superconducting Phsise 

Consider first the effect of the 2DEG on the proper- 
ties of the superconducting phase. In the absence of 
the 2DEG, the superconducting state sustains a "plas- 
mon" excitation associated with oscillations of the phase 
of the superconducting order parameter. This mode is 
gapless in two dimensions with energy ~ \fk at low mo- 
menta fc. This is in contrast to the linear dispersion of 
the sound mode obtained in neutral supcrfluid systems 
(i.e with short-ranged interactions). Upon coupling the 
charged superconductor to the 2DEG, should we expect 
to recover the linear dispersing mode due to the screen- 
ing of the long-ranged part of the Coulomb interaction? 
As shown below, the answer is no. On the time scale 
relevant to the plasma oscillations, the slow diffusive mo- 
tion of the electrons in the 2DEG is unable to screen 
out the Coulomb interaction. Thus the plasmon survives 
essentially unmodified. 

To establish this result, note that in the superconduct- 
ing state, the phase of the Cooper pairs can be treated as 
a classical variable. We can then write down the classical 
equations of motion for the phase (/)(r, i) which obeys the 
Josephson relations: 

dt<^ir,t)=2eU{r,t) (1) 
-2e/9.V<^(r, t) = Jc (2) 

where Jg is the supercurrent, U (r, t) is the potential and 
Ps is the superfluid stiffness (assumed constant). From 
charge conservation in the JJA layer, we have {pc is the 
Cooper pair density): 

dtpc + V-Jc = (3) 

The potential U (r, t) is gencratcxl by both the Cooper 
pair charge density as well as the electronic charge den- 
sity in the 2DEG (pe): 
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U{r,t) = Jd\'V{\r-r'\){p,{r',t)+pe{r',t)) (4) 

(5) 

where the interaction at large distances falls off as Fdr- 
'^'l) ~ \r-r'\ convenience, we have neglected the 
small but finite vertical separation between the JJA and 
the 2DEG. 

Finally, to completely describe the dynamics of this 
system, we need the equations of motions that govern 
charge motion in the 2DEG: 



Je{r,t) = -<Je\7U{r,t) - PVpe 

= dtPe + V • Je 



(6) 
(7) 



where V is the diffusion constant for electrons, and Je 
is the current density in the 2DEG. From equations 
(1,3,4.6), we obtain the equation of motion for the phase 
variable, which is most conveniently written in terms of 
its Fourier components (j){q,uj) defined by: 



d^rdte-'''-''e''^*(t>{r,t). 



The equation of motion is: 

. W2 1 



(8) 



(9) 



where the effect of the 2DEG is simply to replace the 
Coulomb interaction by the effective interaction V^ff. 
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1 



(10) 



Veff{k,uj) V{k) -ioj + m'^ 
The Fourier transform of the 1 /r Coulomb interaction in 

27r 
|fc| 



two dimensions V{k) takes the form V{k) = Stt in the 



long wavelength limit. 

In the absence of the 2DEG {a,, = 0), we obtain the 
usual two-dimensional plasmon mode: 

ojpiik) = sJ{2efpsV{k)k'^ = V8^^\/fc (11) 

Notice that with the 2DEG present, in the static limit 
a; = 0, the Coulomb interaction is screened: 



v;//(fc,c. = o) = i/(|l + |) 



(12) 



For situations that are not static, we cannot automati- 
cally assume this screened Coulomb form for the interac- 
tion, rather, the dynamics of the 2DEG will have to be 
taken into consideration. 

The plasma mode is obtained by solving: 



uj(kf 



1 



(2e)2 V,ff{kMk)) 



Psk'^ = 



u}{k) = u}pi{k) — i{'KCFe)k 



(13) 



Now, since the oscillation frequency of this mode goes as 

the square root of the wave vector, it is much larger (in 
the long wavelength limit) than the damping rate which 
depends linearly on wave vector . Thus, as claimed above, 
the plasmon is a well defined excitation and retains its 
Vk dispersion despite the presence of the 2DEG. 

This problem of the phase dynamics deep inside the 
superconducting state in the presence of a metallic gate, 
was examined earlier in [11] where, in contast to the 
present work, the plasmon was found to be overdamped 
at long wavelengths. We believe this result to be in er- 
ror; in particular, the form of the dielectric function used 
there (equation (20a) of reference [11]) is valid only in 
three dimensions and not applicable to the two dimen- 
sional geometry of the current problem. 



B. Nature of the Phase Transition 

We now consider the nature of the zero temperature 

phase transition between the superconducting and insu- 
lating states of the JJA. In this section we present some 
general arguments that strongly suggest the presence of 
a fixed line of critical points controlling the generic (in- 
cluding disorder) superconductor-insulator transition of 
charged particles in the presence of the metallic plane. 
Note that in the presence of quenched disorder the su- 
perconductor to insulator transition will be a continuous 
transition in two dimensions [15]. 

Consider a system of bosons in two dimensions inter- 
acting via a 1/r Coulomb potential and coupled to a dif- 
fusive 2DEG in the manner described. Such a system 
can be described in terms of a charged bosonic field min- 
imally coupled to the time component of a gauge field 
[Aq). The spatial components of this gauge field can be 
ignored since we are in the extreme type II limit, at zero 
external field. Here, we will examine the action governing 
the gauge field, which takes the following form in imagi- 
nary time (a detailed derivation is in the next section): 



Sa = IJ d-\d-ioV;j}{k,uj)\Ao{k,u;)\'' 

2-^J^'kdru;^-^\Ao{k,ojr 

^ kdroj— — —\Ao{k,uj)\^ 



2e* 



U + Dk'^' 



(14) 



(15) 



which, in the long wavelength {k small) limit reduces to: 



where e* = 2e is the charge of the Cooper pairs, and we 
have introduced the compact notation d^u = 

The first term on the right hand side of equation (15) 
gives rise to the Coulomb interaction while the second 
term arises from the coupling to the diffusive 2DEG. No- 
tice that these terms contain non-analytic functions of 
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k,uj. If wc pcrfom a Rcnormalization Group transforma- 
tion on this system, by integrating out the boson and 
gauge field modes that lie within an energy-momentum 
shell of finite thickness, then additional terms of this form 
cannot be generated. Hence the scaling of these operators 
arises solely from rescaling of space-time and of the fields. 
Since the rescaling of the gauge field is fixed by the re- 
quirement of gauge invariance, the scaling dimensions of 
these operators can be easily derived. Gauge invariance 
requires that the gauge field must scale as the inverse 
time, to preserve the minimal coupling form of the ac- 
tion. Hence its Fourier transform in two dimensions, Aq, 
scales as: 



[Ao{k,w)] 



(16) 



where L has dimensions of length. Thus the Coulomb 
term has dimensions: 



(fkduj\k\\Ai){k,uj)\^] = l}- 



(17) 



where z is the dynamical critical exponent ([time] = L^). 
Thus, the coupling strength of the Coulomb interaction 
e*^ has the renormalization group flow: 



de* 



dln{h) 



\z-l) 



(18) 



when modes between [A, A/6] are integrated out. Thus, 
a fixed point that has a finite value of the Coulomb cou- 
pling necessarily also has z = 1 [7] . This can also be seen 
from the following simple argument. The 1/r Coulomb 
interaction in two dimensions, is non-analytic in Fourier 
space. Therefore, the inverse linear power law of the 
interaction cannot be renormalized. In addition, if the 
fixed point has a finite value of the Coulomb charge, then 
the energy must scale as 1/r which immediately implies 
that the dynamical critical exponent is unity. 

We can now apply a similar analysis to the second term 
in equation (15). In the long wavelength limit, the Vk'^ 
term in the denominator is irrelevant, so long as the dy- 
namical critical exponent for the theory satisfies z < 2. 
This condition will need to be checked at the end of the 
calculation to ensure consistency. Using the known scal- 
ing dimension of the gauge field, we find the interesting 
result: 

[j d'kdu:^AAo{k,Lj)\'']=JJ' (19) 
J M 

Not only is this operator marginal from its scaling di- 
mension, it is always strictly marginal! This is an exact 
statement, arising from gauge invariance and the non- 
analytic form of the gauge field action, and is not just 
the engineering (tree level) dimension of this operator. 

The presence of a strictly marginal operator in the 
fixed point action will in general lead to a line of fixed 



points. This conclusion should obtain if the renormal- 
ized Coulomb charge remains non-zero at the transition. 
Universal critical properties (such as exponents) will vary 
continuously along this line, and hence will depend on 
CTg. As argued below, the Coulomb charge is expected 
to be non-zero for the generic superconductor-insulator 
transition that occurs in the presence of disorder, from 
a superconducting phase to a gapless Bose-glass phase. 
This points to the existence of a fixed line in the presence 
of the 2DEG. 

In the absence of the metallic gate, the argument for 
a finite value of the Coulomb charge at the superconduc- 
tor to Bose-glass transition was reported in [7]. There, 
the singularity in the long wavelength behaviour of the 
compressibility (arising from the Coulomb interaction) 
on both sides of the transition was used to argue that 
the Coulomb charge remains finite at the transition it- 
self. Although this argument fails in the presence of the 
metal due to screening in the static limit, a different ar- 
gument can be made in this case to obtain information 
about the generic transition. Imagine that at this transi- 
tion, the Coulomb charge does not take on a finite value, 
but instead flows to zero. In this limit the gauge fleld 
will decouple from the bosons, which then only experi- 
ence short-ranged interactions. The transition will there- 
fore be in the dirty boson universality class for which it 
is known that z = 2 [16]. However, this scenario would 
not be consistent with the flow equation (18) which in- 
dicates that near this fixed point the Coulomb charge 
grows and renders it unstable. Hence we are led to con- 
clude that the Coulomb charge will remain non-zero at 
the transition between the Bose-glass and the supercon- 
ductor. This leads to two possible scenarios. In the first 
scenario, the Coiilomb charge takes on a finite value at 
the transition, as in the case without the metal. Then, a 
line of fixed points all with z = \ will obtain. In the sec- 
ond scenario, the Coulomb charge flows away to infinity, 
and the gauge fluctuations arc completely controlled by 
the second term on the right hand side of equation (15). 
Such flxed points coifld be stable for z > 1, and once 
again we are lead to expect a line of fixed points, but 
with z > 1. For consistency, we will also require z < 2 to 
justify our assumption in deriving (19) where the term 
containing the diffusivity which is quadratic in the mo- 
mentum was neglected in comparison to the frequency. 
The question of z > 2 fixed points in the presence of the 
metallic 2DEG will not be discussed in the present work, 
and is left for future study. 

Further support for a line of fixed points is provided 
by our calculations on a model of the superconductor- 
insulator transition with particle-hole symmetric disor- 
der, where these scenarios are explicitly realised. The 
critical properties of this model are explored within a 
double epsilon expansion, and both a line of fixed points 
with z = 1, as well as a separate line of fixed points 
with z > 1 are found, corresponding to the two scenarios 
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described above. 

The situation in the clean problem, where the transi- 
tion is between the superconductor and the Mott insu- 
lator is, however, different. For this special case, despite 
the presence of the strictly marginal operator, none of 
the universal quantities are affected, since the Coulomb 
charge flows to zero at the transition, which is hence con- 
trolled by a single (3D XY) critical point. We will never- 
theless consider this problem first, since it will allow us 
to erect our formalism in a technically more simple case. 
The properties of this model are then considered within 
the large- A'' and e expansions. Subsequently, we consider 
a modified model, that includes particle-hole symmetric 
disorder, which behaves more like the generic case. This 
model is studied using the e,er expansion [12] [13] which 
reveals that indeed the universal properties of the tran- 
sition depend on the value of the 2DEG conductivity CTe, 
and hence a fixed line results. 



III. SUPERCONDUCTOR-INSULATOR 
TRANSITION IN THE CLEAN LIMIT 

In this section we describe a model of an array of 
Josephson junctions coupled to a diffusive 2DEG via 
Coulomb interactions. Wc ignore the presence of disorder 
and assume that the juncitions are at integer Cooper-pair 
filling. (Note that in the clean limit, at non-integer filling, 
the insulating phase would generally break translational 
symmetry - a direct transition from the superconductor 
to such an insulator is expected to be first order). 

Consider an array of superconducting islands at posi- 
tions labelled by r. We model each of these islands as an 
0(2) quantum rotor; this is a legitimate approximation 
when there are on average many Cooper pairs present on 
each island. The deviation of the boson charge from the 
background value is denoted by n^. Then, the number 
of Cooper pairs on each site ^ is conjugate to the phase 
on each island cbr . 



e* 



(20) 



Coupling the charge fluctuations in the JJA layer to the 
diffusive 2DEG via Coulomb interactions as described 
earlier, we obtain: 

H = —Ej cos{(l)r — (pr') 

+ 1 + ^ne{rW{r - r')(n^' + Sne{r')) 



2DEG 



(21) 



where Srig is the fluctuation of electronic charge density 
in the 2DEG, the Coulomb interaction at large distances 
is taken to be V{r — r') ^ l/27r|r — r'\ (the factor of 
2-K is introduced for later convenience). i?2D£;G 



Hamiltonian of the diffusive 2DEG, without the Coulomb 
interaction between electrons. We now recast this as a 
path intergal over the bosonic {n,4>) and fermionic (V'e) 
degrees of freedom. It is also convenient to introduce a 
Hubbard- Stratonovich decoupling of the long-range in- 
teraction using an auxilliary field Aq. (As will be clear, 
Aq may be interpreted as the scalar component of a gauge 
field.) This leads to the following representation for the 
finite temperature (= 1//3) partition function: 



Z = I VcpVnVAoVtpee' 



(22) 



where the phase satisfies the boundary conditions 
(j>r{P) = <l^r{0) + l-KYTir with any integer, and 



S = Sb + Sint + Sa[Ao] + S2DEGbl'e 



(23) 



-Ej ^ cos(<^,(r)-0,'(r))] (24) 

<rr'> 

Sint = \ Y^Ao{r,T){nr{T) + 5ne{r,T)) (25) 
= ^ / Y.Mr,r)V-\r - T')A,{r' ,t). (26) 

^ rr' 

Here S2UEG ^^'^ action correspoding to the diffu- 
sive 2DEG in the absence of the Coulomb interaction, 
= dr and implies the matrix inverse of the 
potential. Following the pioneering work of Hertz [17] on 
quantum phase transitions in fermionic systems, we inte- 
grate out the diffusive electrons in the 2DEG to derive an 
effective action for the bosons. The contribution arising 
from this proceedure is: 

Z[Ao]= lvi,ee-^°--oe-i^l^r^°(^^^^^+'^(^'^^^ (27) 

= 2[0]<e^l^^ ^°("+'"^ >2DEG (28) 

Performing a cumulant expansion for the average, and 
taking into account < (5n(r, r) >2deg= we have: 

- ^ / < Snr{T)6nr'{T') >so^^^ Ao{r,T)Ao{r',T') 

■^6 Jt,t' ^ 

+ ... (29) 

As in Ref. [17], we shall neglect the three and higher body 
interaction terms that are generated to obtain: 

ZXA,\ - e- nO(.-.',.-.')Ao(.r)^o(.',r') ^3Q^ 

where 
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n°(r-r',T-r') = ^2 < ^"r(T)^n^'(T') >2DEG ■ (31) ^^^^ bounary conditions <?!.^(/3) = </.^(0) + 27rm^, 
6* where is integer. In the superconducting state, we 

We thus obtain: can replace the cosine in the action (38) by a quadratic 

gradient term. The quadratic action for the phase fluc- 
Z = j 'D(j)'Dn'DAo'De~^''-'f tuations in the superconductor is: 



^eff = f '^r{dr'l> + ^o) ' EjY, COs{<i>r - K) (32) 

+^ Mr,r)V-f){r -t',t- r')^o(r',r') (33) 



5«e = y (T'fc^ra; [^y-}(fc,a;) + {Ejo?)k''\\4>k^\ (39) 



where 



y,7}(r,r)=F-i(r)5(r)+nO(r,T). 



(34) 



where a is the microscopic cutoff, and the integral over 
modes in momentum space is restricted to |fc| < \/^n ja. 
This combined with the 'Lindemann' condition [A is a 
For a diffusive metal the Fourier transform of 11° is given constant of 0(1)): 
by: 



n°(A;,a;) = 



(35) 



where <7e is the conductivity and V is the diffusivity of ^^^^ fo^^i^g -^^^-^^ equation for the critical 

the metal. Thus, the Fourier transform of the effective 



potential is: 



1 



1 



fc2 



Veffik,u) V{k) 

= \k\+a,_ 



|w| +£>fc2 
fc2 



Iwl + 



(36) 
(37) 



point Ey. 



A= dr kdrui {—:r[\k\ + ae 



U+Vk"^ 



]+Eja^k^)- 



Note that in the static (w = 0) limit, the Coulomb in- 
teraction is screened as can be seen in equation (12). 
However, the question of screening is more subtle when 
it comes to dynamical phenomena (as we have seen in 
the case of the plasmon). 



A. Shift in Phase Boundary 

In this section we study the dependence of the zero 
temperature superconductor-insulator phase transition 
point on the 2DEG conductivity. Let us call the critical 
Josephson coupling Ej; for Ej > Ej the system is in the 
superconducting phase. We would like to find the depen- 
dence of Ej on the 2DEG parameters. A rough estimate 
of this phase boundary is obtained by starting out in the 
superconducting phase and applying an analogue of the 
Lindemann criterion; i.e. we ask for what value of the 
parameters are the superconductor phase fluctuations of 
0(1). 

It is convenient for this analysis to rewrite the partition 
function for the sytem in terms of the phase variable: 

2^ J P^exp -i J ^\<Tio^[V-/fik,io)]\4>ko.\'' 



+ Ej 



COs{<f)r — I 



(38) 



(40) 



Inspecting the integral it is clear that the supercon- 
ductor is stabilized in the presence of the 2DEG, i.e. the 
phase transition point E'j{ij,-^ is reduced compared to its 
value in the absence of the 2DEG: E'j{(T,) < E'jiO). This 
is schematically depicted in Figure (2). Qualitatively, 
this effect is to be expected, since the screening arising 
from the 2DEG contributes to weakening the Coulomb 
interaction, and hence enhances the stability of the super- 
conduting state. Also, it has been argued [18] that cou- 
pling the phase fluctuations to external degrees of free- 
dom ('dissipation') leads to reduced phase fluctuations 
which would also stabilize the superconducting phase. 
This effect has been observed in the experiments of Rim- 
berg et al. [4] who were able to cross the insulator to 
superconductor phase boundary just by increasing the 
conductivity of the 2DEG. Experiments on Josephson 
junction arrays with chromium shunt resistors (that do 
not go superconducting), have been reported in [20]. A 
phenomenological model of that system can be written 
down which, in the superconducting phase, would be ex- 
actly the same action as in equation (39) above without 
the Vk^ term, with the conductance being the con- 
ductance of the chromium resistor. In that experiment 
too, it was found that increasing the value of increases 
the stability of the superconducting region. 
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FIG. 2. The Superconductor-Insulator phase boundary in 
the presence of a 2DEG (conductivity ae). The 2DEG stabi- 
lizes the superconductor, as described in the text. Thus, the 
critical value of the intergrain Josephson coupling Ej at which 
superconductivity is established, decreases with increasing Ce- 



tions - the large-A'^ technique and the e expansion. As 
shown below, in both these approximation schemes, the 
Coulomb interaction is found to be marginally irrelevant, 
even in the presence of the metal. As a result, the gauge 
field decouples from the bosonic degrees of freedom at the 
critical point, and the presence of the strictly marginal 
operator (the coupling to the 2DEG) does not affect the 
critical properties of the system. The transition is there- 
fore controlled by the same fixed point as in the clean 
case with short-range interaction and without any cou- 
pling to the metal. The agreement between the large- 
N and e expansions is strong evidence that neither the 
Coulomb interaction or the coupling to the 2DEG af- 
fects the critical properties in the clean commensurate 
superconductor-insulator transition. 

In the generic case, away from this idealized limit, we 
however expect that the combination of the Coulomb in- 
teraction and the coupling to the 2DEG will lead to a 
fixed line controlling the transition. This is illustrated in 
the next section by introducing disorder into the model. 

1. Large N Calculation: 



B. Critical Properties: Clean Case 

In order to study the critical properties of this system, 

it is convenient to reformulate the action in equation (38) , 
by rewriting it as a "soft spin" model, where, after coarse 
graining, the rotor "spin" e~"^ is replaced by a complex 
field ?/'• This automatically takes into account the non- 
trivial boundary conditions on the phase </>(/?). The finite 
temperature partition function can then be written as a 
path integral over these fields: 



S^= j dV(iTi|(5,-ao)V|' + c|VVr+r|V;|2 + ||V>|^ 



2e 



1 f 2 

-j^ / ^ kdrui [\k\ + Ue 



||Ao(fc,c.)|2 (41) 



Since we are interested in the critical properties of this 
model, we can safely ignore the Pfc^ part in the denomi- 
nator of the second term in (41). Also, we will scale out 
the bare sound wave velocity c that appears in the ac- 
tion above by a suitable rescaling of the imaginary time. 
All velocities will therefore be measured in units of this 
bare velocity. For convenience, we continue to use the 
same symbols for the scaled quantities, but these factors 
will have to be restored when comparing with the exper- 
iments. 

In this section we consider the effect of the 2DEG 
and the Coulomb interaction on the properties of the 
superconductor-insulator transition in the absence of dis- 
order. This is studied within two different approxima- 



The effect of the 2DEG on the critical properties of 
the clean problem can be considered within the large- 
approximation. We remain vn D — 2 spatial dimensions, 
but generalize the model of a single species of boson to 
one that has flavours of complex boson fields ■;/',: 

N 

S^= (frdTj2\{dr - iAo)^if + IVV-il' + ro|Vi|^ + 

i=i 

ElV'il']' (42) 



u 



Sa 



N 
2^ 



(t^k(tu[\k\+ae-r-;]\Ao{k,uj)\' 
\u>\ 



(43) 



As argued previously, the entire renormalization of the 
Coulomb charge e* comes from rescaling the fields and 
the momenta. This leads to the following exact equation: 



de 



*2 



dlnb 



= {z- l)e 



*2 



(44) 



while the term containing cTg is exactly marginal. Here, 
we will derive an expression for the dynamical exponent 
(z) at the critical point - within a large- A^ approximation. 

In the N ^ oo limit, the (saddle point) solution to 
the above problem has = 0, and the critical point, at 
this order, is identical to that of the short-ranged model. 
We therefore need to go to next order in 1 /N to see the 
effect of the gauge field coupling . In calculating z to 
order 1/N at the critical point, we only need to consider 
the self energy diagram shown in figure (5a) but within 
this approximation the gauge propagator Ga is replaced 
by its RPA form: 
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*-2 



[|fc| +cre 



(45) 



as shown in Figure 5c. Therefore we find the dynamical 
exponent is given by: 



l = ^/(e*^ae) 



N 



(46) 



.77/2 

Jo 



de 



sin2 6'(l-3sin2 



27r2 e*-2(l + CTe tan6») + %^ 



(Details of this derivation arc relegated to Appendix A). 
In the absence of the metal and for small Coulomb charge 



we find: /(e* ^ 0, CTe = 0) = - 



o*2 



327r' 



and implies 



that a weak Coulomb interaction is marginally irrele- 
vant at the short-ranged fixed point in the absence of 
the metal. Even in the presence of the metal, it is found 
that /(e*^,(Te) < 0, and hence the Coulomb interaction 
continues to be marginally irrelevent in the large- ap- 
proximation. Thus, the only fixed point solution to the 
flow equation (44) is e* = 0. At this flxed point, the 
gauge field coupling may be ignored and the presence 
of the metal does not afii'ect any of the universal critical 
properties. Therefore in the large- A approximation de- 
spite the presence of the metal, the transition is in the 
universality class of the short-ranged model. 



on integrating out high energy modes. Thus equation 
(44) again obtains for the fiow of e* while the term con- 
taining fje is exactly marginal. Here, we will calculate z 
to first order in e, the relevant diagram is shown in Figure 
7a. We find (see Appendix B 1 for details): 



o*2 



71" Jo 



sin2 0(4sin2^-l) 

do ; 

1 -h (Te tan 



(50) 
(51) 



where A{ae) is normalized so that ^4(0) = 1. In that 
limit, CTe — 0, we recover the flow equation of Ref. [21], 
so that the Coulomb interaction is marginally irrelevant 
since the right hand side of (44) is negative. Even for 
non zero Ce, it can be shown that A{ae) > and thus 
the Coulomb interaction is marginally irrelevant despite 
the presence of the metal. The fixed point value of the 
Coulomb interaction then is e* = 0. As a result, the 
gauge field is decoupled from the bosons at the critical 
point, and the presence of the strictly marginal operator 
does not affect the universal critical properties. 



IV. LINE OF FIXED POINTS IN A MODEL 
WITH DISORDER 



2. The e Expansion for the Clean Problem 

In the dimensionality of interest D = 2, the quartic 
term [i/'l* is relevant by power counting. One way to 
control the RG flows is via the standard technique of 
considering the problem close to its upper critical dimen- 
sion, the dimension at which the quartic coupling is also 
marginal by power counting (_D = 3). Therefore, we write 
the effective action (41) in = 3 — e dimensions as: 

S = S^ + Sa (47) 
S^ = yrf^rdr|(at-ao)^|' + |VV'|2+r[^|2 + ||V'|^ (48) 

Sa = ^ J ct''Mu{\k\''-' +aej^^)\Ao{k,Lu)\^ (49) 

Notice that the action for the gauge field is written in 
a general dimension D, so as to always produce the 1/r 
form of the Coulomb interaction (if o-g = 0). The term 
containing is continued to general dimensions so that 
it remains an exactly marginal operator. 

This problem, in the absence of the metal (ue = 0) 
was considered in [21] where the Coulomb interaction was 
found to be marginally irrelevant. Here too we will inves- 
tigate the flow of the Coulomb coupling, but in the pres- 
ence of the metal. For D < 3, the non-analj^ticity of the 
Coulomb term in momentum space once again ensures 
that it is not renormalized by terms that are generated 



In this section, we present evidence to support our 
claim that the critical properties for the generic transition 
(as opposed to the idealized clean commensurate limit in 
Section HI above) are controlled by a line of fixed points. 
We study a model of the disordered superconductor- 
insulator transition in the presence of Coulomb interac- 
tions and a metallic plane. This is in contrast to the 
zero disorder transition considered above. (Note that 
even in the absence of the metal the clean commensurate 
superconductor-insulator transition with Coulomb inter- 
actions has behaviour different from the generic case). 
We consider instead a slightly more realistic model that 
contains a special form of disorder that preserves the par- 
ticle hole symmetry of the problem. In the absence of the 
metal, such a model was recently studied by Herbut [13] 
who found a. z = 1 fixed point with a finite value of 
the Coulomb charge. Here we add the coupling to the 
metal, and first demonstrate that a line of fixed points 
with a finite value of the Coulomb charge and z = 1 can 
be obtained. Subsequently, a line of fixed points with 
2 > 2; > 1, for which the Coulomb charge flows to infin- 
ity, is also obtained. 

The model we consider is essentially the one in equa- 
tion (47) but with a coupling to disorder V, included: 

S = S^ + Sa + Sdis (52) 
S^ = J d^rdT(|(at-ao)V|' + |VV|2+r|^|2 + J|^|4) 
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Sa 



1 



-D, J- /I, |_D-1 , _ M J , ,m2 



Cf fc(fw(|fc|'' 



2e*2 

= / d^rdrV(r)|V(r,r)|2 



ae-^)\Mk,u;)\' 

lUJl 



Notice, that the coupling to disorder preserves the par- 
ticle hole symmetry of the model (the action is invari- 
ant under the particle-hole transformation: -0 — > -0*, 
Aq — > —Aq). For quenched disorder, the random poten- 
tial V is independent of imaginary time, and therefore has 
a big effect on the physics of the system. In the following, 
it is convenient to generalize it to a random variable that 
is correlated in dimensions; in other words we assume 
the following statistics for the random potential: 

< V > = (53) 

< v(7^)v(7^') > = ws"+^-'- {r - r') (54) 

where the label R runs over both space and time com- 
ponents. For quenched disorder e,- = 1. However, it will 
be convenient to study the model in the limit of small 
Ct which will allow one to control the effects of disorder. 
This approximation was introduced in Ref. [12] to study 
the problem of bosons in the presence of disorder and 
short-ranged interactions. 

In order to perform the disorder average, we use the 
standard replica technique. We introduce n copies of the 
fields ^ — > "00, a = 1 . . . 71, average over different realiza- 
tions of the disorder potential, and finally take the ti — > 
limit. Averaging over disorder generates the following 
term which has a non-local interaction in e-r dimensions 
between the replica fields. 

n 

S = J2SBm + SA + S'^,, (55) 

a=l 

SL-~Y ^ Jd^+^Rd^+^R'S'^+'-^^iR-R') 



a, 13=1 • 
|2| 



|0„(i?)|^|^^(i?')|^ 



(56) 



We now consider performing an RG transformation on 
the above action, by integrating out the short wavelength 
modes 1/A < a < 6/A, and then rescaling momenta, fre- 
quency and fields to obtain the flow equations for the 
various couplings. The calculation is controlled by work- 
ing in Z? = 3 — e space dimensions, and as mentioned, 
with disorder that is correlated in dimensions. We as- 
sume that £,67- are small and derive the flow equations to 
lowest order in these quantities. Details of this derivation 
are relegated to Appendix B. After a redcflnition of the 
coupling constants: 



q = 
A = 

w = 



U 

8^ 
W 
2^ 



the flow equations take the following form: 

^ = (--1)? =(y- 3^)9 (58) 
d\ B , 5 o C , 11 

dh^ = ^^+2^'^ -2^ -T'^ 

^^^[e + er)W + l;W'^ + ^q^W -2\W (60) 
dlno 8 2 

and the coupling constant i = ^ is strictly marginal 
and does not flow. The parameters A,B,C that appear 
in the flow equations are functions of Ue and have been 
normalized so that ^(0) = B{Q) = C(0) = 1. They are 
given by the following expression: 



A{a,) 
where 



2 r^'^ 

- ddsm^diAsin:^ 0-1) f{ae,0) (61) 

Jo 



IT 12 



dO sin^ Of{ae,d) 



- / desm^ef{ae,0) 

Jo 



1 



1 + (Te tan ( 



(62) 
(63) 

(64) 




(57) 



FIG. 3. The results of the e, Et expansion showing the fixed 
lines obtained in the presence of the marginal coupling 1/g 
which is proportional to the 2DEG conductivity. Only the 
fixed point value of the Coulomb charge e* is plotted. The 
coupling constants e* and g are measured in units of e, where 
er — e has been assumed. The z = 1 fixed points are obtained 
for < , while z > 1 fixed points are obtained for all 
g. The fixed points with both e* and g large are not to be 
trusted since gauge field fiuctuations are big in that case, and 
are shown with dashed lines. The limiting cases shown as C 
and SR, correspond to the Coulomb fixed point found in [13] 
and the short-ranged fixed point found in [12] 
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A. Fixed line with z = 1 



3.5 



We now solve for the fixed points of tiie flow equations 
above, assuming that the Coulomb charge e* takes on a 
finite value at the transition. The coupling contants at 
the fixed point are proportional to e but depend both on 
the ratio e^/e and on cTg. In the following wc will as- 
sume Cr/e = 1, a choice that is consistent with our aim 
to study the = e = 1 limit. We evaluate the functions 
A,B,C numerically for various values of CTg and find the 
fixed point values of the couplings which are quantities of 
order e that depend on the value of cTg. It turns out that 
fixed point solutions of this kind exist for CTc in the range 
< ae < 1.99. While it is convenient to solve for the 
fixed points in terms of (Te, they are more appropriately 
labelled by the value of the marginal coupling g, which 
can be obtained from the relation: ^ = ^^in . When 
expressed in terms of this variable, the fixed points exist 
in the range | < ^ ~ (572e)^^. In fact, as shown in fig- 
ure 3, for every value of g in this range, there are a pair 
of fixed points. The fixed points with a smaller value of 
the Coulomb charge, evolve in the limit ^ — * (which 
corresponds to removing the metallic plane) to the fixed 
point studied by Hcrbut in [13]. 

All of the fixed points obtained are found to be stable 
within a linear stability analysis; although the eigevalues 
are found to be complex, they have negative real parts. 
Thus we have a stable line of fixed points for a range 
of the 2DEG conductivity. The dynamical critical expo- 
nent z for all of these fixed points is unity. Other critical 
exponents can be expressed in terms of the fixed point 
couplings: 



= eql{ue)B{ae) 



while the exponent v can be expressed as: 



(65) 



2.5 



1.5 



0.5 



1 n"^ ^ 

FIG. 4. The exponents r] and v for the z = \ fixed line 
(lower branch in Figure 3) controlling the S-I phase transi- 
tion. Along the horizontal axis is plotted the marginal cou- 
pling 1/g which is proportional to the 2DEG conductivity. 
For this plot we have set e = er = 1. From the dependence 
of the exponents on the conductivity it is clear that the criti- 
cal properties are significantly affected by the presence of the 
2DEG. 

The remaining exponents can all be derived from a 
knowledge of these, if hyperscaling is assumed. The 
Chayes inequality v >2/D is satisfied by all these fixed 
points. A plot of these two critical exponents as a func- 
tion of the marginal coupling 1 / g (which is proportional 
to the 2DEG conductivity) is shown in Figure (4). Only 
the exponents for the lower branch of the z = \ fixed line 
in figure 3 are shown. The exponent v is found to increase 
with increasing 2DEG conductivity. If this qualitative ef- 
fect that is present in the e expansion persists unchanged 
in the physical problem, then this would predict that the 
critical region for a finite sample would increase as the 
2DEG conductivity increases, which could be experimen- 
tally accessible. We return to this point in the following 
section. 

This calculation therefore, provides us with an explicit 
demonstration of the general claims of this paper - that in 
the presence of the Coulomb interaction and a metal, the 
dirty superconductor to insulator transition can be con- 
trolled controlled by a line of fixed points, parametrized 
by the metal conductivity. 
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-9c(<^e) 



) (66) 



B. Fixed Line With z>l 



We now turn to the question of whether there are any 
where the function D{ae) (which has been normalized to stable fixed points in this model with e* = oo. For clarity 
be unity in the absence of the 2DEG, D{0) = 1) is given we rewrite the gauge action as: 
by: 

Sa = ^J cr^fc<ra;|fc|^-i|io(fc,c^)r 

D{(Te) = - I ' desm^26f{(Te,e) (67) +^ f dt^kdtuj^\AQ{k,uj)\'^ (68) 

71" Jo ^9 J M 
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where again g has bccns used to denote the coupling 
strength of the marginal term, and is related to the bare 
couplings by 5 = e*'^ /(jg. Under the Renormalization 
Group transformation, g docs not flow. Now, if the renor- 
malized Coulomb charge were to flow to inflnite strength, 
then the gauge fluctuations would entirely be controlled 
by the second term in the gauge action in equation (68). 
In order that this is a stable flxed point we need to en- 
sure that .2 > 1, so that /3(e*) > at the flxed point. 
In this limit, the flow equations can be obtained from 
the expressions (58 - 64) by formally taking the limit e* , 
cTe — > oo but holding the ratio e*^/cre finite and equal to 
g. Defining: 

9 

yields the following flow equations: 
d\ , 1 , 5 , , 1 



dln5 
dW 
dlnb 



TT Z 



while 



(e + er)W + Iw"^ + -gW - 2\W 

8 TT 



8 37r 



(69) 
(70) 

(71) 



Assuming that g is a constant of order e, and working 
with e = Ct-, we can solve for the fixed points in this the- 
ory to first order in epsilon, as a function of the marginal 
coupling g. Fixed points are found to exist for all such 
values of g, and they are all found to have z > 1. A lin- 
ear stability analysis of these fixed points demonstrates 
that they are stable (eigenvalues of the stability matrix 
are found to be complex, but with negative real parts). 
Thus, a stable fixed line of critical points with z > 1 can 
be obtained within this model. 

In order to derive this model we have assumed that the 
dynamic critcal exponent satisfies z < 2. This allows us 
to retain just the frequency part in the denominator in 
passing from equation (41) to equation (49), which is the 
form used in all the subsequent analysis. In the epsilon 
expansion, the dynamical critical exponent differs from 
unity by an amount of order e. Since e is considered the 
smallest quantity in the problem, the condition on the 
dynamical critical exponent z < 2 is formally satisfied 
within our calculation. 

These results are summarized in Figure 3, where only 
the couplings e* and g for the fixed points is plotted. The 
z> 1 fixed points are found for all values oil/g (which is 
proportional to the conductivity of the 2DEG). However, 
for very small values of 1/g they are not to be trusted, 
since we assume g to be 0(e). Therefore, these fixed 
points are shown with dashed lines. In the limit of very 
large 2DEG conductivity {1/g large), the fixed point ob- 
tained is the same as for the problem with short-ranged 
interactions and no metallic plane, which was considered 
in [12] (this is made explicit in the Appendix B4). 



V. QUANTUM HALL TRANSITIONS WITH A 
METALLIC PLANE 



In this section we consider the effect of a metallic plane 

on the quantum Hall transition, which, along with the 
superconductor to insulator transition, has been one of 
the most intensively studied quantum critical phenom- 
ena [24]. Most theoretical studies of the quantum Hall 
transition between integer quantum Hall states have used 
a non-interacting model of electrons in a random poten- 
tial subject to a magnetic field. While it is believed that 
short-ranged interactions between the electrons are irrel- 
evant [6], the same is not true of the long-ranged Coulomb 
interaction which is present in the physical problem of 
interest. It has been suggested that the non- interacting 
electron universality class for these quantum Hall tran- 
sitions may be partially recovered by screening of the 
long-range Coulomb interactions by a nearby metallic 
gate. Increasing the screening in this way was shown 
experimentally to induce a change of variable-range hop- 
ping transport, away from the transition, from the Efros- 
Shklovksii behavior expected with Coulomb interactions 
to the Mott behavior expected for noninteracting elec- 
trons. [26] So one conjecture would be that a nearby 
metallic gate simply changes the transition to the short- 
range universality class. 

However, as we have seen in the case of the 
superconductor-insulator transition, the effect of the 
metallic gate may well be more involved, and it cannot 
automatically be assumed that its sole effect is to screen 
the Coulomb interactions. Indeed, the general arguments 
in Section II B (relying as they did only on gauge invari- 
ance and the non-analytic form of the gauge field action) 
strongly suggest that a similar fixed line may also control 
the generic quantum Hall transition in the presence of a 
metallic plane. 

To provide some calculational evidence for this sug- 
gestion, we consider a simplified model used to describe 
the quantum Hall transition with interacting particles, 
although the effects of disorder, which are believed to 
be crucial for the physical transition, are ignored. The 
model describing the critical point consists of a mass- 
less Dirac fermion coupled to a Chern-Simons gauge field 
[27] and interacting via the Coulomb force [14]. This 
theory can be used to model both integer and fractional 
Hall transitions by changing a statistics parameter. The 
Euclidean action for this model, in the absence of the 
metallic gate, takes the form: 

S = Sf -\- Sgauge 

f ^ - i 

Sf = d^xdr ^ ^/'m[7o(9o - ^^"'o) 

m=l ^ 
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(72) S„ = l^%^.Mk,.)e,kM-k,-u.) 



S„ = j ^\^.[MKo.)e.,kM-k,-u:) 



+ -Y Vcoui (k) \eijkiaj{k,uj)f] (73) 

where the 7^ are the Dirac matrices in three dimen- 
sional EucUdean spacetime (one representation in terms 
of the famihar Pauli matrices is [70 = icrs] jj = icrj]), 
e* is the Coulomb charge, 9/N the statistical angle and 
Vcoui{k) = 1/k corresponds to a Coulomb potential of 
In order to exert some control over the physics of 
the critical point, N species of fermions were introduced. 
When N = 1 the model describes charged particles at the 
critical point between an insulator and a quantized Hall 
state with Hall conductivity axy = = 

this is a model of the integer transition, while 6 = 2it 
describes a Mott insulator to supcrfluid transition. In 
reference [14] the model was considered in the limit of 
large N. The results of those investigations to order 1/A'' 
were as follows. The critical properties are controlled 
by a line of fixed points parametrized by (the statis- 
tics is not renormalized, as expected). For 6 < 1.24, the 
Coulomb interaction is found to be marginally irrelevant, 
while for > 2, it flows to strong coupling. In between, 
1.24 < ^ < 2, the critical points are found to have a finite 
Coulomb charge and z = 1. 

We now consider a system undergoing an insulator- 
quantized Hall state transition in the presence of a metal- 
lic 2DEG. The geometry is taken to be identical to figure 
1, except that in place of the J J A, a layer undergoing the 
quantum Hall transition is present. The metallic 2DEG 
below it is assumed to have sufficiently high density of 
electrons that we can ignore the effect of the applied mag- 
netic field, and treat it as a diffusive metal. 

To model the quantum Hall transition, we use the 
anyon Mott transition model described above, but now 
the particles in the quantum Hall layer interact with the 
electrons in the metal, via the Coulomb force. Once again 
we integrate out the diffusive electrons in the metal to ob- 
tain an effective action for the quantum Hall transition. 
Again, only the part of the response that is quadratic 
in the gauge fields is retained, while higher order terms 
are neglected. It is convenient for our purposes to write 
down this action in a somewhat different form from that 
shown above. An additional temporal gauge field Aq is 
introduced, and the Euclidean action can then be written 



as: 



S — ^Dirac ~t~ ^gauge 

N 



Soirac = / d^xdr ^ V'm[7o(9o /^("^o + ^0)) 

■' m=l ■V-'^ 



where 



V^f){k,i^) = \k\+ae 



(75) 



(76) 



Note that gauge invariancc and the non-analytic form 
of the terms entering the gauge part of the action (75) 
ensure that the statistical angle 6 as well as the term 
containing the conductivity (7^ are not renormalized, and 
the Coulomb charge flows according to the equation 
(44). Thus, there are two strictly marginal operators 
parametrized by 9 and (Te/e*^- We already know that 
the statistics parameter generates a fixed line of critical 
points. We would now like to address the question: under 
what circumstances, at a fixed value of 6, will the 2DEG 
also give rise to a fixed line? We will address this ques- 
tion within the large- A/' approximation. As we have seen 
in earlier cases, if the Coulomb interaction turns out to 
be marginally irrelevant, then there is a single fixed point 
controlling the transition (at a fixed 6) despite the fact 
that the term containing the 2DEG conductivity cr,-. is 
a strictly marginal operator. This is because the gauge 
field ^0 decouples from the problem in this limit. We 
would therefore like to calculate the flow of the Coulomb 
charge to 0(1/A^), in the presence of the metal. This can 
be done following reference [14], but with the modified 
action (75). For example, for the integer transition = 0, 
the flow equation for the Coulomb charge in the presence 
of the metal can be easily worked out {w = e*^/16): 



dw 
dlnb 



Jo 



sm 



1 -|- w sin^ + ae tan 6 



(77) 



+7j(5i - -^ai)]V'm 



(74) 



the right-hand side clearly is always negative and hence 
the Coulomb interaction at this integer transition is 
marginally irrelevant even in the presence of the metal. 

The flow equations for general 6 and de are compli- 
cated and we do not attempt to write them out here. 
Instead we will argue that at least for small enough (t,, 
there is a range of 6 flxed points for which, if wc fix 6, 
a fixed line parametrized by (Te can be obtained. First 
consider the case of (Te = where a finite value of the 
Coulomb charge was found for fixed points in a certain 
range: 1.24 < 6 < 2 [14]. Imagine fixing a particular 
value of 9 inside this range. Then, for sufflciently small 
(Te > 0, a fixed point with finite Coulomb charge is still to 
be expected. The strictly marginal operator contained in 
the CTe term will affect the universal properties of the fixed 
point by inducing a fixed line at this particular value of 
6. Since the renormalized Coulomb charge is expected to 
remain finite along this fixed line, z = 1 will obtain. In 
addition to these fixed points, we can consider the effect 
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of the metal (with (Je small) on the fixed points for 9 > 2. 
There, the e* ^ oo behaviour in the model without the 
metallic plane, is expected, in the presence of a metallic 
plane, to lead to a line of z > 1 fixed points. These fixed 
points should also be accessible within the large N tech- 
nique; although the Coulomb charge flows to infinity the 
gauge field fiuctuations are small and controlled by the 
marginal coupling e*^/(7eiV. 

Thus, in this model of quantum Hall transitions, as in 
the superconductor-insulator transitions discussed previ- 
ously, the proximity of the 2DEG can give rise to a line of 
fixed points controlling the critical properties and hence 
affect the phase transition in a non-trivial fashion. 



VI. COMMENTS ON LOCAL DISSIPATION 
MODELS 

In this section, we briefly comment on previous the- 
oretical studies of Josephson arrays with "local dissipa- 
tive baths" whose results bear superficial similarities with 
those in the present paper. It is important to first realize 
that we have focused on a physical situation where the 
origin of the dissipation is very clear. It is due to the cou- 
pling between the Cooper pairs and the gapless diffusive 
electrons in the metal. In contrast, Ref. [9] postulated 
the presence of dissipation due to some unspecified local 
degrees of freedom. That paper also suggested a fixed 
line controlling the transition between the superconduct- 
ing and insulating phases, once the dissipation exceeds a 
critical strength ao — 2/3. The dissipative term for each 
superconducting grain is assumed to be 

S, = „/*..-<Mp)!. (78) 

This model is not expected to describe the experiments 
in Refs. [4,5]. The zero diffusion constant limit of our 
model resembles superficially the local dissipation models 
considered in [18], which in turn differs from (78) in that 
the local dissipation is coupled to phase differences across 
junctions. Moreover, in our case the crucial boundary 
condition on the phase in the imaginary time direction 
(38) can differ from that derived for more phenomeno- 
logical models of local dissipation [18,19]. Furthermore, 
the presence of long ranged interactions is crucial to the 
physics discussed in this paper. 

Finally we point out a problem with the approxima- 
tion used in [9] to derive the presence of a fixed line in 
the absence of both long-ranged interactions and disor- 
der. We show in Appendix C, that when exactly the 
same approximation is applied to the 3D anisotropic XY 
model, one erraneously concludes the existence of a fixed 
line in that problem, in which the critical properties are 
actually controlled by a single fixed point. 



VII. IMPLICATIONS FOR EXPERIMENTS 

In this section we will consider the physical conse- 
quences of the field theoretical results derived earlier. 
We first address a quantitative question - at what value 
of the 2DEG conductivity are the critical properties of 
the superconductor-insulator transition significantly af- 
fected? For the z = 1 fixed points, this will clearly oc- 
cur when the value of the marginal coupling (l/g) is of 
the same order or larger than the renormalized Coulomb 
coupling (1/e*^). Reinstating the physical parameters, 
this yields the following condition on the 2DEG con- 
ductivity in SI units: af''^* > {^)2csKeo, where k€o 
is the permittivity constant and Cg is the bare speed of 
sound in the superfluid (in the absence of Coulomb in- 
teractions). If we make the naive assumption that the 
renormalized Coulomb force is roughly the same as the 
bare Coulomb force e* w e*, then we can obtain an es- 
timate for the 2DEG conductivity at which a significant 
effect on the critical properties is felt. The speed of sound 
is given by (cg ~ ^^/EjEc), where a is the distance be- 
tween superconducting grains and Ej and Ec are the 
Josephson coupling energy and the charging energy of 
a grain. Typically, at the transition Ej ~ Ec (if the 
2DEG conductivity is not too large) and so we obtain 
Cg ^ Eca/h. Finally, substituting the form of the charg- 
ing energy for a grain of radius R, Ec = e*'^ /A-Kne^R, we 
obtain that the nondimensional conductivity should be 
2cKeo ~ ^{a/R). This is roughly of order the universal 
conductivity, if the grain size is approximately the same 
as the distance between grains. Hence, we expect the ef- 
fect of the 2DEG to be substantial once its conductivity 
is of order (6.4kn)~^. This regime is readily accessible 
in experiments. 

An alternate estimate of this conductivity for the z = 1 
fixed line can be derived from the epsilon expansion solu- 
tion to the model with particle-holc symmetric disorder. 
There, we know that the exponents of the z = 1 fixed 
points are strongly affected once the marginal coupling 
is of the order of l/g,; = 1.7 x 10^'^e. Setting e = 1 and 
using the relation crf''^" = (27r)^i, we obtain a^^^'' = 

0.01^ ~ (600fcf2)~^. This corresponds to a smaller 
value of the conductivity than the earlier estimate, be- 
cause the Coulomb charge at the fixed point is strongly 
renormalized. For the z > 1 fixed points, the marginal 
coupling controls the gauge field fluctuations. However, 
once the 2DEG conductivity greatly exceeds the quantum 
of counductance l/g » 1 =^ Ue > (^ti")^ = (I.IA;^)"^ 
the critical exponents approach those of the short-ranged 
model without a metallic plane. 

We now briefly consider how the results for the flxed 
line of critical points obtained in this paper may be tested 
experimentally, both for the superconductor-insulator 
transition as well as the quantum Hall transition. Ide- 
ally, one would like to investigate the universal critical 
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properties at the transition, such as the critical expo- 
nents and the conductivity at the transition, as a func- 
tion of the 2DEG (metaUic ground plane) conductivity. 
A fixed hne of critical points implies that these properties 
are a function of where the phase boundary is crossed. 
Although these various 'universal' quantities would vary 
along the fixed line, there would be universal relations 
amongst them, since they can all be expressed as a func- 
tion of a single parameter, the strength of the marginal 
coupling. We also note the caveat, that for large val- 
ues of the 2DEG conductivity, there will be region about 
the critical point where the system will behave like the 
short-ranged interaction model, before ultimately cross- 
ing over to the true fixed point behaviour very close to 
the transition. 

Finally, we note that in a finite sized system (such as 
the JJA in reference [4] which is a 40x40 system) the 
transition is rounded off when the correlation length ex- 
ceeds the linear size of the system. The width of this 
region (Ag), in terms of a control parameter g that tunes 
the transition, is given by the exponent f {Ag ~ L~^l^ ^ 
where L is the linear system size). The double cpsilon ex- 
pansion described earlier, on a model with particle-hole 
symmetric disorder that exhibits the superconductor to 
insulator transition, revealed a strong variation of v along 
the fixed line. Such a variation could presumably be de- 
tected by studying the width of the transition in these 
finite sized systems. 



VIII. CONCLUSIONS 

In this paper, we have studied the effects of a proxi- 
mate metallic plane on various two dimensional localiza- 
tion transitions with an emphasis on the superconductor- 
insulator transition. Such a metallic plane has two prin- 
cipal effects — it provides a mechanism for screening the 
long-ranged Coulomb interaction, and it is a source of 
dissipation due to the gapless diffusive electrons. The in- 
terplay of these two effects leads to interesting physical 
phenomena. Perhaps the most interesting result in this 
paper is the possibility of a fixed line with variable critical 
exponents controlling the transition. In addition, right at 
the transition, a temperature independent conductivity 
that will vary along the fixed line is expected. 

Our results are of direct relevance to experiments prob- 
ing the superconductor-insulator transition in Josephson 
junction arrays in the presence of a metallic plane [4] . So 
far these experiments have not probed the universal scal- 
ing regime near the transition. We hope that our results 
will focus future experimental work on this regime. 

A similar scenario may be expected to hold for quan- 
tum Hall transitions in the presence of a proximate metal- 
lic plane. A line of fixed points would again imply vari- 
able critical exponents and a temperature independent 



value for the diagonal conductivity at the transition. Re- 
cently, quantum Hall systems with nearby metallic planes 
have been prepared [26] and it would be interesting to ex- 
perimentally investigate the effect of the metallic plane 
on the critical properties at the quantum Hall transitions. 

One limitation of our field theoretic approach is that, 
following Hertz [17], wc have integrated out the gapless 
diffusive electrons to obtain an effective action for the 
critical bosonic Cooper pair degrees of freedom. As with 
other problems involving quantum phase transitions in 
fermionic systems, a more satisfying theoretical approach 
would keep both gapless fermionic and bosonic modes as 
part of the effective theory and treat them on equal foot- 
ing. Unfortunately, such an approach is not available at 
present. 

Another assumption that has been made which is 
needed to justify an approximation used through this 
paper is that the fixed points satisfy z < 2. We leave 
it to future work to investigate whether superconductor 
to insulator transitions in the presence of Coulomb in- 
teractions and a metal can be controlled by fixed points 
with dynamical exponents 2; > 2, and if so what their 
description is. 

We conclude by noting that our results show that 
even a simple Coulomb interaction between the Cooper 
pairs and the gapless fermionic degrees of freedom 
has a profound effect on the universality class of the 
superconductor-insulator transition. Should the fixed 
line found in this paper survive the inclusion of processes 
where the Cooper pair can decay into two electrons, it 
will lead to a temperature independent though appar- 
ently non-universal conductivity at the superconductor- 
metal transition. 
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APPENDIX A: LARGE N CALCULATION OF 
Z - 1 IN THE CLEAN CASE. 

In this section we describe in some detail how the dy- 
namical critical exponent z can be calculated within the 
large- iV approximation, in the absence of disorder. 
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We consider the problem directly in D = 2, but with 
species of bosons coupled to the scalar component of 
a gauge field as described by equation (43). The bare 
propagators for the gauge field and the bosons are: 



0*2 



N \uj 



(Al) 
(A2) 



(a) 



(b) 



= \/\/\/^ 



(c) 



FIG. 5. To evaluate the 0(1/N) correction to the dynami- 
cal critical exponent; only diagram (a) needs to be considered, 
where the solid line is the boson propagator at criticality and 
the dashed-wavy line is the RPA form of the gauge propaga- 
tor, defined in (c), that appears to this order in the large- iV 
approximation. 

Since it is only the dynamics and coupling of the gauge 
field that distinguishes space and time, any contribution 
to ^; — 1 must come entirely from this interaction. We 
have already noted in Section III B 1 that to leading or- 
der in 1 /N (saddle point) the gauge field decouples from 
the problem, and the dynamical critical exponent stays 
equal to unity at this stage. Therefore, we calculate this 
quantity to 0(1/N). First, we recognize that the gauge 
field propagator will be dressed by boson density fluctu- 
ations as shown in Figure 5a, so that we obtain the RPA 
form: 



(A3) 



where II" is the contribution from the pair of boson bub- 
ble diagrams in Figure 5c 



n°(fc,c^) = J 



ct Kctngo{K,Q.) 



[{uj + 2nfgo{k + K,uj + n)-i] (A4) 

= - , " (A5) 
16 Vfc^ + 

In order to calculate {z—1) to 0(1/N), we need to evalu- 
ate the self energy contribution arising from the diagrams 
in Figure 5b. We will be interested in isolating terms in 
the self energy that are of the form: 



E{k,uj)^r]k^ ln(\/A:2 + uj^/A) + rj^to^ ln(\/fc2 + uj^/A) 

(A6) 



where A is the ultraviolet cutoff. These contributions will 
modify the boson propagator since: 



g-\k,io) = [g°{k,u;)]-'^ - E(fc,a;) 



(A7) 



and thus affect the power laws that characterise the boson 
correlators at criticality. The dynamical critical exponent 
is then given by: 



z-l = ^{Vw-V) 



(A8) 



Notice, that since we are only interested in contributions 

of the form (A6), which vanish in the zero frequency limit, 
we only need to calculate S(fc, uj) — S(0, 0), therefore the 
second diagram in Figure 5b, which gives a momentum 
independent contribution, can be ignored. Evaluating 
the resulting diagram: 

E{k,u) - S(0,0) = 
j dr'^Kdn [{2u + nfg^{k + K,w + fl) 

-n'^g'^{K,n)]gAiK,n) (A9) 

we can obtain the correction to the dynamical critical 
exponent to 0(1/N), whic is: 

z-l = l/(e*2,ae) 



/(e*^C7e) 



d9 sin2 0(l - 3sin2 6') 
2^e*-2(l + aetan0) + 2^ 



Note, that this can be shown to imply that z — 1 > for 
all e* > 0, (Te > and hence the only fixed point is at 
e* = 0. 



APPENDIX B: THE e, EXPANSION 

We consider the problem in £) = 3 — e spatial dimen- 
sions, with particle-hole symmetric disorder correlated 
in e-r dimensions, given by the action in equation (55). 
We consider integrating out all modes with wavevectors 
in the range [A/6, A]. The additional terms generated 
lead to a renormalization of the terms in the action; 
the quadratic derivative terms are now multiplied by the 
scale factors Zk and the other terms acquire fac- 
tors Zr, Zu and Z^r- Note, in dimensions D < 3, the 
Coulomb term, as well as the term arising from coupling 
to the 2DEG in the gauge propagator are unchanged on 
integrating out high frequency modes since they are non- 
analytic in frequency wave-vector space. On rescaling 
the momenta k' = bk to restore the cutoff to^A and 
also rescaling uj' = Vu) and the fields ip' = Z^'^ and 
A'q = b~^AQ (the gauge field scaling follows from gauge 
invariance) to obtain a new set of couplings at this scale: 



16 



w' = b'+^-'+'-z^^z^w 

*'2 



r' = b^Z-^Zrr 



(Bl) 
(B2) 
(B3) 
(B4) 



In the following, we will look for fixed point solutions to 
the above flow equations, assuming that e, Ct are small. 
In that limit, wc arc justified in looking for fixed points 
where the couplings are small, and hence a perturbative 
evaluation of the flows is feasible. 



The Feynman rules derived from the action (55) are 
shown in Figure 6, where the boson propagator at criti- 
cality is given by = {k^ + ui^) and the gauge propa- 
gator is given by Q^^ = e*~^(fc^~^ + (Tej^)- 




(a) 



q,(B 



(b) 



1. Calculation of z and 77; Flow of e*. 

Here we calculate the change in the quadratic deriva- 
tive terms {Zk, Z^^) on integrating out the high momen- 
tum modes. In practice, this is most easily done by in- 
tegrating all modes below a cutoff A, and extracting the 
log divergent contributions to the boson propagator Q: 



g-^i{k, = + _ ^^2 in(Vfc2 + coy A) 



(B5) 



More details on how these two proceedures are related to 
each other may be found in [23] . Here, we just note that if 
integrating all modes below the cutoff A gives rise to the 
additional term in equation (B5), then the contribution 
on just integrating modes in a shell betwen [A/6, A], can 
be found by subtracting the same term but with A ^ A/6 
which yields: 

{i]k'^ + rj^uj'^) In 6 

. These can then be related to the quantities we want to 
calculate: 




Zk = l-|-?7ln6 
Zoj = 1 + riujlnb 

q+Q, 
2a)+i2 



^AA/ 



(B6) 
(B7) 




FIG. 6. Feynman rules used in the e, er expansion. The 
vertices arising from the jV'j* interaction, the gauge inter- 
action and disorder generated interaction are shown. The 
dashed disorder line only transfers the D + l — er spatial com- 
ponents of the momentum. 



(c) 

FIG. 7. The diagrams (a) and (b) contribute to renormal- 

izing the the quadratic derivative terms in the action, while 
all three diagrams contribute to the renormalization of the 
quadratic term r. 

We are now in a position to calculate the quantities 
in (B7). Since we are interested in the momentum de- 
pendent contribution to the self energy, wc only need to 
calculate w) — E(0, 0). The two diagrams of interest 
then are shown in Figure 7a and 7b. We work to lowest 
order in e and therefore the integrals will be performed 
in D = 3, and the terms that depend on the logarithm 
of the cutoff will be extracted. The contribution from 
gauge fluctations (Figure 7a) is: 

s(fc,w) - E(o,o) = J (t^KttngAiK^n) 



[{2lo + nfg{k + K,w + n)- n^g{K, n)] (bs) 



which yields: 



n d6'4-9cos^l 
Jo IT I +(Te tan 6 
e*2 p ^6* (4 cos- (9 - l)coti-0 
677^ Jo TT l + (Te tan 9 



"^0 + 4003^0 



(B9) 
(BIO) 



Similarly, the contribution arising from the disorder di- 
agram in Figure 7b can be calculated. Since the disorder 
induced interaction is independent of spatial momentum 
transfer, T,{q, 0) = T,{0, 0), the disorder contribution to r] 
vanishes. However, at finite frequency transfer we have: 



S(0,a;)-S(0,0) = -a;2w j 



(Bll) 



which, to lowest order in Cr can be evaluated to give 

W 



w 



8772 



(B12) 
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putting this all together we have 



from this the flow equation for the Coulomb charge is 
easily derived: 



*2 



(B13) 
(B14) 



after redifining the couplings = l::^ and W — — 



we get equation (58). 



2. Flow Equations for u: 

We begin by calculating the renormalization of the 
quartic term (u) arising from the integration of high 
frequency modes (Zuu). Diagrams Figure 8a-g con- 
tribute, and we can write the result as ZuU = M-|-(Fig. 
8a+..+g)ln6. We will evaluate these diagrams and ex- 
tract the log divergent parts; the external lines are as- 
sumed to carry zero momentum. 




(a) 



(b) 



(c) 





(e) 




(f) 



(g) 



FIG. 8. The diagrams contributing to the renormalization 
of the quartic interaction. 

The contribution from the first diagram can be found 
in any standard reference on the e expansion for short- 
ranged model [23] , and is given by: 

5 «2 
^^^-28^^ 

Next, contributions arising from the gauge interaction 
(Figure 8b-e) can be evaluated. Details about the com- 
binatorial factors that appear for these diagrams are not 
repeated here as they as they are discussed in [22]; the 



only difference for our calculations is that the gauge prop- 
agator is modified due to coupling to the 2DEG. The 
result for these diagrams then is 



8b-t-8c-t-8d ^ -^[- 



g*4 
2^2 1 



tt/2 



de- 



(l + (Te tan 61)2^ 
while the graph in Figure 8e contributes: 

„*2 ^ /•7r/2 



8e- 



ue 



47r2 



de- 



cos 



. tan^ 



The diagram in Figure 8f which involves the disorder can 
be computed in the following way. There are ^2 dia- 
grams of the type shown, and hence the contribution is: 



8f ='^C2uW 



J 



3uW 
1^ 



Finally, the diagram in Figure 8g makes no contribution. 
Since the external lines are at zero frequency, and the 
disorder interacton, being independent of time, does not 
change the frequency, the coupling to the gauge field van- 
ishes. 

Putting together these results, as well as the factors 

arising from rcscaling the fields (Bl, B6, B14) and re- 
defining the couplings as in (57), we obtain the flow equa- 
tion (59). 



3. Flow Equations for Disorder Coupling 

We flrst consider the renormalization of the disor- 
der term (W) arising from the integration of high fre- 
quency modes (Zy^W). Diagrams Figure 9a-d con- 
tribute, and we can write the result as ZyyW ■ 



V^+(Fig. 



9a+..-|-d)ln6. 




(b) 




(c) 



(d) 
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FIG. 9. The diagrams contributing to the renormahzation 
of the interaction induced by disorder. 

To evaluate the diagram in Figure 9a we note there 
are two identical processes that contribute, with the bo- 
son lines parallel or anti-parallel. Similarly, there arc two 
equivalent processes associated with the diagram in Fig- 
ure 9b, with the disorder induced interaction decorating 
cither the left or right boson line. These diagrams all 
make the same contribution; putting them together: 



9a-F9b = 



^4 

a p 



27r2 



Similarly, for the diagram in Figure 9c there are two 
equivalent contributions, since the gauge interaction can 
decorate either the left or the right boson line. Hence, 

9c = 2w I dr^Kdrnn^gA{K,n)g^{K,n) 



'■W A 



47r2 



Tr/2 



d6- 



cos^ 6 



1 + (Te tan ( 



Finally, the the diagram in Figure 9d also has two equiv- 
alent realizations, and the boson self interaction loop can 
be inserted in four ways, which yields: 



9d = m{-u/4) 



^4 

a p 



-W 



47r2 



Putting together these results, as well as the factors aris- 
ing from rescaling the fields (B2, B6, B14) and redefining 
the couplings as in (57) , we obtain the flow equation for 
the disorder coupling (60). This completes the derivation 
of the flow equations (58 - 60). 



4. The z>l Fixed Line 

The fixed points in the presence of a metallic plane in 
the case e* ^ oo are obtained from the equations (69, 
70) assuming that the marginal coupling g is 0(e). This 
yields the fixed point values of the couplings, (in terms 
oix{g) = \{^g + Ae + lle,)): 



Ac 

Wr. 



7 _ 

-g'^+x'^{g)) 



(B15) 
(B16) 



the couplings are 0(e) and the dynamical critical expo- 
nent is given by: 



1 + 



9_ 

37r 



(B17) 



These solutions exist for all 0(e) values of g, and have 
z > 1, which ensures that e* = oo does not flow. In 



the linit g ^ (which corresponds to a large value of 
the 2DEG conductivity) we recover the short-range in- 
teracting boson case, in the absence of a metallic plane 
(studied in [12]). Note, the flow equations (69, 70 with 
5 = 0) however differ from those in [12]. This difference 
is a result of treating the time direction as distinct in 
the presence of Coulomb interactions [13], which leads us 
to scale the measure [divdh^] as L~^^^ rather than as 
as in [12] where only the e,- dimensions are 
assumed to scale with the dynamical exponent. There- 
fore when we continue the fixed point obtained in the 
presence of Coulomb interactions to the ^ = fixed point, 
it differs in details from the result obtained in [12] but 
agrees with the short-ranged fixed point quoted in [13]. 



5. The Critical Exponent v 

To calculate the exponent we need the flow equa- 
tions for the coefficient of the quadratic term (r). In 
fact, we only need those terms in the flow equations that 
themselves contain r. The diagrams that then need to be 
considered are those in Figure 7a,b and c, evaluated with 
zero momentum and frequency in the external lines. In 
these, the terms proportional to r that contribute to its 
renormahzation are: 

7a = r j dt^Kdtng'^{K,n)n'^gA{K,n) 



dO- 



cos 



27r2 TT Jo l + cTg tan 6 
7h = rW j dr%g^{p) 
W 

finally, evaluating the diagram in Figure 7c, we note that 
there are four ways to pick boson lines to make a loop 
and hence: 

7c = 4r(- J) / ,f%g\p) 



—r 



87r2 



Combining this with the effect of rescaling the fields, as 

given in equation (B4) and rescaling the couplings ac- 
cording to (57), we obtain the following flow equation for 
this relevant coupling. 



dr 
d\nh 



= [2 - A + 



W D 
T 12 



q^]r + . .. (terms independent of r) 

(B18) 



Where the coefficient _D is a function of Ce and is defined 
in equation (67). At the fixed point this leads to the 
following equation for the critical exponent v: 
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W, D 2 



(B19) 



APPENDIX C: THE ANISOTROPIC XY MODEL 



since these fixed point couplings are all small [0(e)], we 
can invert this expression to obtain equation (66). In 
the limit cTe = 0, we can compare this answer with ear- 
lier work. The coefficient of the term differs from the 
expression in [13], but agrees with that in [22] 
The Exponent v for the z > 1 Fixed Line 
For this fixed line, u can be derived from equation (66) 
by formally letting g^,fTe oo but holding the ratio 
q^/ae = g finite. This yields: 



, = l + i(A,(^)_A_^(^)) 
2 4^ cvy; 3^ ^ \yjj 



(B20) 



6. Stability of the Fixed Points: 

We would like to make sure that besides the relevant 
variable r, the fixed point is stable to perturbing away in 
any of the other couplings. A linear stability analysis can 
be performed on the fixed point at (g^, VFc, Ac) by move 
slightly away (g^ + 6q^, Wc + 6W, Ac + i5A) and asking 
how the deviations evolve under resclaing. This leads to 
the equation, 



d 



dlnb 



Sq^ 

6W 
6X 



where 



d„.p{q^) 

d,2 3iW) 



dwPiq') 
dwl3{W) 




dx(3{q^) 



L 



(B21) 



(B22) 



where the derivatives of the beta functions are evaluated 
at the fixed point. If the eigenvalues of Mc are all neg- 
ative, then the perturbation will die out on running the 
RG, and the fixed point is stable. 
Stability of z = I Fixed, Points 

Wc have numerically evaluated these eigenvalues along 
the z — 1 fixed line parametrized by < (Tc < 1.99, and 
find that all these critical points are stable (we have as- 
sumed € = St which is consistent with the limit of e = 1, 
Er = 1 that we are finally interested in). A pair of the 
eigenvalues, it turns out, is complex (and conjugate to 
each other), but they have negative real parts. The least 
negative eigenvalue ranges from ~ — 0.5e (at CTc = 0) to 
0.3e (near = 1-99). 

Stability of z > 1 Fixed Points Here, the stability ma- 
trix we need to consider is just a 2 x 2 matrix since the 
condition z > 1 which is satisfied by these fixed points 
ensures 1 /e* = remains stable. The eigenvalues of this 
matrix were evaluated numerically (assuming e = e,-) for 
the full range of g (assumed 0(e)) and were found to be 
complex but with negative real parts. Thus, all these 
points on the z > 1 fixed line are found to be stable. 



This appendix shows that the arguments in 
fixed line for the dissipation term 



-.J 



drdr 



,(0( T)-0(rO)^ 

1-^ _ ^'|2 



for a 



(CI) 



which unlike in our model do not depend on either long- 
ranged interactions or disorder, predict incorrect behav- 
ior in a simpler case and hence are unreliable compared 
to controlled approximations such as the e-expansion. 

The approach of [9] predicts a continuous set of transi- 
tions for a simple related problem, an anisotropic XY 
model, where only one superconducting transition is 
present. This suggests that phase-squared dissipation 
likely introduces at most one new universality class for 
the superconducting transition. The subtle problem in 
the derivation of the effective action in [9] is that only 
the leading (two-spin) term is kept in the Hubbard- 
Stratonovich transformation used to decouple the junc- 
tions. 

Consider a classical 3D system made of layered 2D XY 
models with interplane coupling cos(0j — ^i+i), where 
i is a layer index. The action after taking the continuum 
limit in the planes is 



SxY 



^ J d^x [J, 



,(V0i)2- J^cos(0i-0i+i)]. (C2) 



Here the couplings are defined to include 13 = 1/ kT. If 
Jz = 0, then the system can have algebraic in-plane or- 
der for Jxy above the Kosterlitz-Thoulcss transition. In 
this algebraically ordered phase, the spin-spin correlation 
{Sr-i ■ S'ra) oc [ti — r2|~", for somc value < a < 1/4 de- 
pending on the in-plane coupling. Introducing a complex 
boson field ip = Sx + iSy to decouple the planes. 



^ = IZ/ rf'a;[J,^(Vei)2 + |Vi-e^^f 

i 

-Jz(V^*V«+i + V'mV'.)/2]. (C3) 

Now the approach of [9] is based on an effective action 
for the full 3D model with nonzero interplane coupling 
Jz, using the above 2D correlations. The quadratic part 
of this action reads in cylindrical coordinates 



/ 



'S'eff = / kdkdcj) dkz {kz 



(C4) 



Now this quadratic part is used in [9] to calculate the 
critical conductivity as support for the claim of a conti- 
nous set of universality classes. 
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FIG. 10. Phase diagram of 3D XY model with direc- 
tion-dependent couphngs Jz^Jx ~ Jy Arrows indicate di- 
rection of projected RG flows. Despite the line of fixed points 
for Jz = 0, all transitions between the long-range ordered (su- 
perconducting) phase and short-ranged (insulating) phase are 
controlled by the single critical point C. 

Even though (C4) reproduces the in-plane two-spin 
correlation function in the absence of interplane coupling, 
it cannot be used at this level of approximation to study 
the transition into the superconductor. The actual phase 
diagram of the anisotropic XY model is shown in Fig. 10. 
There is a single fixed point that controls the supercon- 
ducting transition for all nonzero values of the couplings, 
with unique values of the critical exponents. 

The above suggests that the model considered in [9] 
may not in fact have a continuous set of superconducting 
transitions once a > ao = 2/3. The dissipation generates 
a series of of terms in the ip action; all of these are ir- 
relevant for small a, but the higher-order terms become 
relevant at the same time as the leading term. Weak 
phase-squared dissipation is indeed irrelevant, but it is 
not clear what happens once the dissipation is relevant. 
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